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Abstract 

We show that the maximum fidehty obtained by a p.p.t. distiUation protocol is given by the solution to 
a certain semidefinite program. This gives a number of new lower and upper bounds on p.p.t. distillable 
entanglement (and thus new upper bounds on 2- locally distillable entanglement). In the presence of sym- 
metry, the semidefinite program simplifies considerably, becoming a linear program in the case of isotropic 
and Werner states. Using these techniques, we determine the p.p.t. distillable entanglement of asymmetric 
Werner states and "maximally correlated" states. We conclude with a discussion of possible applications of 
semidefinite programming to quantum codes and 1-local distillation. 

1 Introduction 

One of the central problems of quantum information theory is entanglement distillation ( Q , |^ ) : the production 
of (approximate) maximally entangled states from a collection of non-maximally entangled states. Of particular 
interest are 1-locally distillable entanglement and 2- locally distillable entanglement (the amount of entanglement 
that can be distilled using local operations and a 1-way (2- way) classical channel). Nearly all of the known upper 
bounds on 1- or 2-locally distillable actually apply to a larger class of operations, known as p.p.t. (positive 
partial transpose) operations This motivates our present study of p.p.t. distillable entanglement. 

We study distillable entanglement via a more refined quantity, the "fidelity of distillation" , which measures 
how close one can come to producing a i^-dimensional maximally entangled state from a given input. In 



Theorem 3.1 below, we show that the fidelity of p.p.t. distillation can be expressed as the solution to a certain 



semidefinite program (see for a survey of semidefinite programming). Then any feasible solution to the 



dual problem (Theorem 3.3) gives us an upper bound on fidelity of distillation. 

The rest of the paper is devoted to an exploration of the consequences of this semidefinite program. Section 
4 gives a number of results that hold in general, including a new bound combining the bounds of [lO[ | and and 
a theorem to the effect that maximally entangled states cannot be used to catalyze fidelity of p.p.t. distillation. 
In section 5, we show that the semidefinite program simplifies in the presence of symmetries; in some cases 
(e.g., isotropic states, Werner states), this simplification turns the semidefinite program into a Zmear program. 
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In the case of asymmetric Werner states, this hnear program can be solved exactly, showing that the upper 
bound of 1^] is tight in that case. Section 6 sketches a technique for producing asymptotic lower bounds, which 
we then use to strengthen the hashing lower bound ||^ in the p.p.t. case. We also use this technique to partially 
resolve a conjecture of by determining the p.p.t. distillable entanglement of "maximally correlated" states. 
Finally, in section 7, we consider possible applications of semidefinite programming to the problems of quantum 
codes and 1-local distillation. In particular, using the techniques of section 5, we give a new derivation of the 
linear programming bound for quantum codes , ||] , ||ll| . 



2 Operators, superoperators and operations 

If y is a Hilbert space, we denote by H{V) the space of Hermitian operators on V. We also let 'P{V) C Ti-iV) 
denote the convex cone of positive semi-definite Hermitian operators; we will freely write ^ > i? to mean 
A — B G ViV). A state is then an element of 'P(F) of trace 1. Quantum information theory can be thought of 
as studying the behavior of these concepts under tensor products. 

Given an operator A G Ti.{V (8) W), we define the "partial trace" Try (A) to be the (unique) operator in 
n{W) such that 

Tr(Try(A)B) = Ti{A{B 1)), (2.1) 

for all B G H{W). Similarly, given a choice of basis for W, we can define the partial transpose A^^ by 

Tr(A^»'(B® C)) = Tr(A(B0C*)), (2.2) 

where B G Ti-iV), C G Ti,{W), and C* is the transpose of C with respect to the chosen basis. Both of these 
transformations extend by linearity to non-Hermitian operators as well. 

A positive operator C G ViV W) is said to be "separable" if it can be written in the form 

C^^A, ®B„ (2.3) 

i 

with A^ G P(V), B^ G ■P(VK); in other words, 

C eViV)(g)V{W). (2.4) 
Similarly, C is said to be p.p.t. (positive partial transpose) if 

c er{v (giW)np{v (g)Wf'''; (2.5) 

note that this does not depend on the choice of basis in W. We also recall that every p.p.t. operator is separable: 

v{v) p{w) c v{v (E)W)n v{v ® wf"^ . (2.6) 

A "superoperator" from V to V is a linear transformation from TLiV) to 7^(1/'). The space of superoperators 
can be naturally identified with TL{V ® V')] to a superoperator corresponds the unique operator 51(5') such 
that 

Tr(B*(A)) =Tr(f7(*)(A®B)). (2.7) 
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We also define the adjoint superoperator by 

Tr(A**(B)) = Tr(B*(yl)). (2.8) 

Note that 

= Ti-v{n{W)(A ® 1)) (2.9) 

^*{B) ^Trv'{n{^){l(g>B)) (2.10) 

and, if *i : n{V) H{V') and : H{V') n{V"), then 

r!(*2 o *i) = (*t ® l)(f^(^2)) ^{l® *2)(f^(*i)) - TiV'((f^(«'i) ® ly")(ly ® ^^(*2))). (2.11) 

Of particular interest is the (self-adjoint) superoperator A^-^ A*] in that case, we find 

n{t) = ^{V^ ® Vi)(Vj (g) Vj)^ > 0. (2.12) 

A superoperator is said to be "positive" if '^{A) > whenever A> 0, and "trace-preserving" if ^*(1) = 1; 
equivalently, Try/ (17(4')) = 1. A superoperator is "completely positive" if it satisfies any of the following 
equivalent conditions: 

• (1) ly ® vj/ is positive 

• (2) For all Hilbert spaces W, Iw ® ^ is positive 

• (3) There exist operators Ai G Hom(y, V') such that 

^{A)=Y,AAAl (2.13) 

i 

• (4) For any (some) basis of V, the partial transpose rj(vl>)'"^ is positive semi-definite. 
Clearly 2 => 1, and 3 2 is straightforward. To see 1 =^ 4, it suffices to observe that 

n{^f''^{iv(E)'^)in{t))>o. (2.14) 

Finally, 4 3 follows by taking an eigenvalue decomposition of il(5')^^. Since the operators we will be 

dealing with in the sequel are mostly completely positive, we define f2'(^') = $1(5')'"^', and use this to identify 
the space of superoperators with 7i(V^ Cg) V). Thus the set of completely positive superoperators is identified 
with V{V ® V). An "operation" is defined to be a completely positive, trace-preserving superoperator; we 
denote the (convex) set of operations from V to V by Op{V, F')-Q 

On tensor product spaces, there are several classes of operations of interest, which can be defined in terms 
of the convex sets V and Op as follows: 

^This differs somewhat from the definition of operation given in jgj, in that we are assuming operations to be "non- measuring" , 
but by the main result of that paper, this incurs no loss of generality when studying entanglement distillation. 
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• e-local: Ce = Op{V, V) Op{W, W). 

• 1-local: Ci = {r{V ® V) ® Op(H/, VF')) n Op{V (g>W,V' (g) W). 

• I'-local: Cv = (Op(T^, V^') ® P{W (g, W')) n Op(F (gW,V' (g W). 

• separable: C$ = {r{V ® F') ® r{W (g W')) H Op{V (g>W,V' (g> W). 

• p.p.t.: Cr = Op(F VK, ® M^') n Op{V (g)W,V' (g W'f«««' . 

We also have the class of 2-local operations, defined by allowing arbitrary compositions of 1-local and I'-local 
operations. For a diS'erent approach to defining these classes, see [||. We recall 

C, cCi,Ci. CC2 CC$ cCr, (2.15) 

with all inclusions strict in general. (The class Cg, not discussed in is simply the closure of the class of local 
operations under convex linear combinations (i.e., shared randomness).) 

From a physical perspective, the only natural classes are those of (e, 1, 1', 2)-local operations. The difficulty, 
however, is that in none of these cases do we have an effective way to decide whether a given operation belongs 
to the class; this is especially true in the case of 2-local operations. Thus the class of separable operations is 
important as a simplification of the class of 2-local operations, while the class of p.p.t. operations is important 
as the smallest class containing the 2-local class for which we can effectively decide membership. For instance, 
all of the known upper bounds on 2-locally distillable entanglement are really bounds on p.p.t. distillable 
entanglement; to a large extent this even applies to upper bounds on 1-locally distillable entanglement. Similarly, 
a lower bound on p.p.t. distillable entanglement provides a limit on how far the current methods can take us. 

3 Fidelity of distillation 

For any integer K > 0, we define the "maximally entangled" state ^(K) G Ti.{C^ ® C^) by 

$(if) - -|r!'(lcK) = -i (e.®e,)(e,®ej)^- (3-1) 

l<iJ<K 

Given any other state p, the "fidelity" of p is defined by 

F{p) = Tr(<i>(if)p). (3.2) 

Definition 1. Let p £ P{V W) be a state, and let ii' > be an integer. The "fidelity of ii'-state p.p.t. 
distillation" Fr{p;K) is defined by 

Fr{p;K)^mi,xFi^ip)), (3.3) 

where ^ ranges over all p.p.t. operations from H{V W) to H{C^ C^). 
Remark. We can define F^, Fi, F2, etc., similarly. 
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This is a refinement of the concept of distillable entanglement; indeed, we can define (see ||, |]): 

Definition 2. Let p be as above. The p.p.t. distillable entanglement Dr{p) of p is defined to be the supremum 
of all positive numbers r such that lim„^oo Fr{p®^; [2''"J) = 1. 

Thus a study of -Fr is likely to provide insights into Dr, as we shall indeed find below. 

We first observe that the optimization problem defining Fr can be rewritten as an optimation over operators: 

Theorem 3.1. For any state p and any positive integer K , 

Fr{p; K) = maxTr(Fp), (3.4) 

F 

where F ranges over Hermitian operators such that 

< F < 1, -1/K <F^ < l/K. (3.5) 

Proof. Let ^ be the operation maximizing F(^'(p)) in the definition of Fr(p; K). Clearly, if we compose with 
any operator of the form U ®U , this leaves F{^{p)) unchanged. The same must then be true after averaging 
over U{K) ("twirling" |^). We may thus assume 5* = T o v]/^ where T is the twirling superoperator. We find 

T{A) = Tr(A$(if ))$(i^) + Tr{A{l - <f{K))){l - ^K)); (3.6) 

T{A) must have the form a^{K) + b{l - ^{K)), and since 

Tr(T(A)) = Tr{A) and Tr{T{A)'i>{K)) ^ Tr(A$(i^)), (3.7) 
we can solve for a and b. It follows that 

n{T) = <^>{K) ® $(if) + ^ (1 - <i>iK)) ® (1 - <^>{K)) = n'{T) (3.8) 

But then we compute 

r2'(*) = n'{T o *) (g, i)in'{T)) (3.9) 

= «'*($(ii:)) eg) <^{K) + — i— 5'*(l - $(i^)) eg) (1 - ^{K)). (3.10) 

— 1 

Setting F = ^'*($(i^)), we obtain: 

^F(S $(if) + —^—{1 - F) (1 - $(if)). (3.11) 

— 1 

This operator is positive if and only if F > and {1 — F) > 0. We also find 

n'i^f ^F^(g <i>{Kf + A l - F^) ® (1 ~ HKf) (3.12) 

iv^ — 1 

= ZTT^^/^ + i^$(i^)^)/2 + - F^) ® (1 - K'^{Kf)/2. (3.13) 

Since (1 ± K^{KY) /2 are orthogonal projections, we find that ri(^')'" is positive if and only if 

-l/K<F^<l/K. (3.14) 
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The theorem foUows by noting 

F{^(p)) = Tr{Fp). (3.15) 

□ 



Definition 3. An operator that satisfies the inequahties ( |3.5| ) wiU be said to be primal feasible for Fy{p;K); 
if it maximizes Tr(Fp), it will be said to be primal optimal. 

We will use this result to define Fr{p;K) for all positive real values of K; for an interpretation, see the 



remark following Corollary 4.3 below 



Theorem 3.2. The function Fy is convex in p and concave in 1/K; that is, for < tt < 1; 

Fr{TTpi + il-TT)p2;K) <nFr{pi;K) + {l-TT)Fr{p2;K) (3.16) 
Fr{p; {KiK2)/{^K2 + (1 - tt)K^)) > 7rFr(p;ifi) + (1 - 7:)Frip; K2). (3.17) 

In particular, Fr is continuous in both variables. 

Proof. Let F be primal optimal for Fr{TTpi + (1 — t^)P2, K). Then 

Frinpi + (1 - tt)p2;K) = Tr{F{npi + (1 - n)p2)) (3.18) 
= 7rTr(F/7i) + (l-7r)Tr(i^P2) (3.19) 
<nFripi;K) + il-n)Fr{p2;K). (3.20) 

Similarly, let Fi and F2 be primal optimal for Fr{p;Ki) and Fr{p;K2) respectively. Then ttFi + (1 — tt)F2 is 
primal feasible for Fr{p; {K1K2) / {ttK2 + (1 — 'k)Ki)), thus giving the second inequality. □ 

The above optimization problem is an instance of what is known as "semi-definite programming" (SDP) 
(p[). That is, it involves the optimization of a linear function subject to the constraint that certain operators 
(depending linearly on the variables) must be positive semidefinite. This has several consequences, including the 
computational one that semi-definite programs can be solved in polynomial time (typically polynomial in the 
dimension, although special structure can greatly reduce this). Another consequence is that there is a notion 
of duahty for SDPs. 

For a Hermitian operator we define the positive part Aj^ and negative part A- to be the unique positive 
operators such that 

A+- A_= A, A+A_=Q. (3.21) 

We also define |^| ^ A+ + A^. 

Tiieorem 3.3. For any state p G 7i{V ® W) and any positive real number K , 

Fr{p;K)= min Tiip - D)+ + ^Ti \D^\. (3.22) 
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Proof. Let F be an operator satisfying the constraints above. Then for any operators A, B, C, we have: 

Ti{Fp) = Tr(A) + l/KTr{B + C) 

- Tr((-p + A + B^ - C^)F) - Tv{A{l - F)) - Ti{B{l/K - F^)) - Ti{C{l/K + F^)). (3.23) 

If A > 0, B > 0, C > 0, and 

A>p^{B- Cf, (3.24) 
then the last four terms are all nonnegative, and we have 

Ti{Fp) <Tr{A) + l/KTi{B + C), (3.25) 

and thus 

Fr{p;K)<Tr{A) + l/KTr{B + C). (3.26) 

In fact, by the theory of duality for SDPs, this inequality can be made tight, to wit: 

Frip; K) = min Tr(A) + l/KTr{B + C), (3.27) 

minimizing over operators satisfying the constraints. Upon adding a variable D with D = {B — C)^ , the 
constraints become 

A,B,C>0, A>p-D, B + C = D^. (3.28) 

We thus find 

Fr(p;K) = min | min Tr(^) + 4 min Tr(B + C) 1 . (3.29) 

D \ A>0 K B,C>0 I 

\A>p-D B-C=D^ I 

But we readily see that 

mm Tr(A)=Tr(p -£>)+, (3.30) 

A>'p-D 

min Tr(B + C) = Tr|£>^|, (3.31) 

B,C>0 \ / I I' \ I 
B-C=D^ 

proving the theorem. □ 

Definition 4. An operator D ^T-LiV ® W) such that 

Fr(p; Jf) = Tr(p -£»)+ + ! TV |£>r I (3.32) 

will be said to be dual optimal for Fr(p; K). 



Thus given any operator Z), we obtain bounds on fidelity of distillation, and conversely any such bound can 
in principle be shown by choosing a suitable operator D. For instance, Theorem |3.2| could also be proved as 
follows: 

Proof. If Di and D2 are dual optimal for i^r(pi; K) and i^r(P2; K), then 

Fr(^pi + (1 - t:)p2\K) < Tr(7rpi + (1 - tt)p2 - nDi + (1 - 11)02)+ + ^ Tr \nD^, + (1 - n)D^\ (3.33) 

A 

< 7r(Tr(pi - Di)+ + Tr + (1 - 7r)(Tr(p2 - ^^2)+ + Tr \DI\) (3.34) 

^7rFr{pi;K) + il-n)Fr{p2;K). (3.35) 

Similarly, if D is dual optimal for Fr{p; {KiK2)/{nK2 + (1 — ■k)Ki)), then 

Frip; iK,K2)/{7TK2 + (1 - 7r)Xi)) = Tr(p - + (tt-^ + (1 - ^)-^) Tr (3.36) 

Al A2 

= 7r(Tr(p -D)+ + -^ Tr + (1 _ ^)(Tr(p - D)+ + ^ Tr |i?i^|) 

Al A2 

(3.37) 

>nF{p;Ki) + {l~7r)Fip;K2). (3.38) 

□ 

4 General results 

Lemma 4.1. For any integer d>l, and any K > 0, Fr{^{d); K) = min(l, c?/i4'). 

Proof. For K > d, take F = d/K, D = $((i). For K < d, take F = $((i), D = 0. □ 

Theorem 4.2. For ott.?/ states pi and p2, and any K,K' > 0, 

Fripi;K')Fr{p2;K/K')<Fr{pi®P2;K)<Fr{pi;K/TT\pl\). (4.1) 

Proof. For the first inequality, let Fi and F2 be primal optimal for Fr(pi; K') and Fy{p2 \ K/K')\ then Fi eg) F2 
is primal feasible for Fy{pi ® P2', K), giving the inequality. 

For the second inequality, let D be dual optimal for Fr{pi;K/ Tr |p2 I)- Then, taking F)' = F* (g) p2 , we have 

Fr(/9i ® P2; K) < Tr((pi ® P2) - {D ® ^2))+ + ^ Tr |(i? ® ^2)"^! (4.2) 

A 

= Tr(pi-F)+ + ^^Tr|Fr| (4.3) 
= Fr(pi;X/Tr|pr|). (4.4) 

□ 

In particular, if Fr(p2; Tr \p\\) — 1, then equality holds in this theorem, taking K' = K/ Tr jpj I- Since this 
is true for ^{d), 



Corollary 4.3. For all integers d, all K > Q, and any state p, 

Fv{p®<^{d)-dK)^Fr{p;K). (4.5) 
Remark. This gives us another way to define Fr(p; K) for general K > 0. For rational K > 0, we can define 

Fr{p;p/q)^Fr{p<i^^q);p), (4.6) 

which is well-defined by the theorem. Since the resulting function is nonincreasing in K , there is a unique way 
to extend it to a left-continuous function of K, which must then agree with our earlier definition. 

Another example is when p2 is p.p.t.; then Tr I = 1- We have: 
Corollary 4.4. For all K > 0, any state p, and any p.p.t. state p' , 

Ft{p®p';K)^Ft{p;K). (4.7) 

Corollary 4.5. For any K > Q and any state p, 

min(l,l/if) < Ft{p\K) < niin(l, Tr |p^|/X). (4.8) 
Proof. By the theorem, we have, writing p — $(1) 05 p. 

min(l,l/if) = Fr($(l);if)Fr(p;l) <Fr(p;X) <Fr($(l);if/Tr|pr|) ^min(l,Tr|pi^|/if). (4.9) 

□ 

Asymptotically, the theorem becomes: 
Corollary 4.6. For any pair of states pi, p2, 

Dvipi) + Dr{p2) < Dripi ® pz) < Cr(pi) + loga Tr (4.10) 

In particular, 

Drip®Hd)) = Dr{p) + \og2{d), (4.11) 

and for any p.p.t. state p' , 

Drip ^p') = Drip). (4-12) 

Remark 1. Subtracting Dr{pi) from the inequality, we obtain the bound Drip) ^ log2Tr|p'"l of [^. (But see 



Theorem 4.13 below.) See also [M, for an independent rederivation. 



Remark 2. For other classes of operations, ( 4.11 ) is known only when Dip) > 
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Since Definition ^ maximizes over all p.p.t. operations, we can obtain relations between different values of p 
and (integral) K by composing with appropriate p.p.t. operations. The next two results extend this. We recall 
from that for a superoperator '5, is defined by 

^"^{p) = ^{p^f. (4.13) 

Note that $!(*'") = , and thus * is p.p.t. if and only if vj/ and are completely positive. 

Theorem 4.7. For any state p, any K > 0, and any trace-preserving superoperator ^ such that both and 
are positive, 

Fr{-^{p);K)<Fr{p;K). (4.14) 

Proof. (First proof) Let F be primal optimal for F('I'(p); K). Then ^'*(F) is primal feasible for F{p; K), so 

Fr{^{p);K) = Tr(F*(p)) = Tr(«'*(F)p) < Fr{p;K). (4.15) 

(Second proof) Let D be dual optimal for Fr(p;K). Then 

Fr{^{p);K) < Tr{^{p) - ^{D))+ + 1/KTi \^{Df\ (4.16) 

< Tt^{{p-D) + ) + l/KTr^^{\D^\) (4.17) 

^Tt{p- D)+ + 1/KTt\D^\ (4.18) 

= Frip;K). (4.19) 

Here we used the facts that for a positive superoperator 5* and an arbitrary operator p, 

Tr(M'(p)+) < Tr(v|/(p+)) and Tr(|vl/(p)|) < Tr(vl/(|p|)). (4.20) 

□ 

Remark. It follows from this that we cannot improve on the p.p.t. fidelity by using trace-preserving superop- 
erators such that both 4" and 4''" are positive. In fact, one can show using the techniques of Section || that 
any such operator that produces isotropic output must in fact be p.p.t. 

Lemma 4.8. For any state p, the function KFr{p; K) is nondecreasing in K, while the function 

{KFr{p;K)-l)/iK -I) (4.21) 

is nonincreasing in K . 

Proof. We first consider KFr{p; K). Writing F' — KF, we have 

KFt{p\K) = maxTr(i^'p), (4.22) 

with F' subject to the constraints 

0<F' <K, -1 < F'^ < 1. (4.23) 
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Since increasing K increases the feasible set, the maximum cannot decrease. Dually, 

KFr{p;K) = mini4:Tr(p -£>)++ Tr |L>^|, (4.24) 

which is nondccrcasing in K for any choice of D. 

For {KFr{p; K) - 1)/{K - 1), we proceed similarly; taking F' = {KF - 1)/{K - 1), we have: 

{KFr{p; K) - 1)/{K - 1) = maxTr(i?^V) (4-25) 

with F' subject to the constraints 

-1/{K - 1) < F' < 1, -2/{K - 1) < F'^ < 0. (4.26) 
These constraints become harder to satisfy as K increases, and thus the maximum cannot increase. Dually, 

{KFr{p; K) - 1)/{K - 1) = min(K Tr(p -£))++ Tr \D^ \ - 1)/{K - 1). (4.27) 



But 



-^(i^ Tr(p - D)+ + Tr \D^\ - 1) = -^(i^ Tr(p - D)+ + Tr(i?r) + 2Tr(Z)r)_ _ i) (4.28) 

= -^(JC TV(p - D)+ - Tr(p - D) + 2 Tv{D^)_) (4.29) 

= Tr(p - D)+ + ;^^(Tr(p - D)_ + 2 Tv{D^)_). (4.30) 

This, of course, is nonincreasing in JC, so we are done. □ 

For integer K, this corresponds to composition by the following p.p.t. operations: 
Lemma 4.9. Let Id{f) denote the isotropic state 

w) := fm + ^(1 - m) (4.31) 

of dimension d and fidelity f. If f <l/d, then for all K > 0, 

Fr{Id{f);K) = l/K. (4.32) 

Otherwise, for < K < d, 

Frihif); K) = 1/K + ^^{1 - 1/K), (4.33) 



and for K >d, 



Fr{Ia{f);K) = ^. (4.34) 
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Proof. For the first claim, take F = 1/K, D = Id{f ), at which point > 0, so Tr \D'^\ = 1. For the second 
claim, take 

F = l/K+^^^^il-l/K) (4.35) 

D = ^^{l + d<i>{d)). (4.36) 

Finally, for the third claim, take 

P^m^ (4.37) 

D = Idif). (4.38) 

In each case, the lower bound coming from F agrees with the upper bound coming from D, and thus both F 
and D are optimal. □ 

Remark. In particular, we have Fr{Idif)', d) — max(l/(i, /); the fidelity of an entangled isotropic state cannot 
be increased by p.p.t. operations. 

It is instructive to translate the relative entropy bounds of , in terms of the dual SDP. We recall the 
definition 

S{p\ \<j) = - Tr(p(log2 p - log2 a)), (4.39) 

and the following result: 

Lemma 4.10. Theorem 2.2] Let p,u E ViV), with p a state. For < e < 1 and n E Z+, define 

7„(e) ~- log2 minTr(a®"P), (4.40) 
n P 

where P ranges over projection operators on 1/®" such that Tt{Pp) > 1 — e. Then 

liminf7„(e) >5(p||a) (4.41) 

n—*oo 

limsup7„(e) < --^S{p\\a). (4.42) 

n — >oo t 6 

Remark. In [Q, this is stated only when ct is a state; scale invariance gives the result in general. Also, if both 
p and a are diagonal, we may restrict P to be diagonal as well; this is just the analogous result of classical 
information theory. 

We then have: 



Theorem 4.11. For any state p and any p.p.t. state a, 

Drip) < S{p\\a). (4.43) 
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Proof. We need to show that for any x > S{p\\a-), 

HmsupFr(/9®";2^") < 1. (4.44) 

n — >oo 

Choose y between x and S'(p||o'), and consider the dual SDP bound with 

D = 2«"(T®". (4.45) 
Then D is p.p.t., so 1/ KTr \D^ \ = 2^y~^^^ 0; the first term is bounded below f by the following lemma. □ 

Lemma 4.12. Let p and a he arbitrary states, and let y be a nonnegative real number. Then 

lim sup Tr(p®" - 2^'"ct®")+ < 1 (4.46) 

n — *oo 

whenever y > S{p\\a). 

Proof. Let Pn{y) be the projection onto the positive part of 

p®"-2^"cr®"; (4.47) 

then we need to show that 

Fn{y) Tr((/7®" - 2^"a«")P„(y)) = Tr(p«"P„(y)) - 2^" Tr(a®"P„(2/)) (4.48) 

is bounded below 1. Fix e, and consider the statement Fn{y) > 1 — e. For this to be true, we must certainly 
have 

Tr(p«"P„(2/)) > 1-e (4.49) 

Tr((7®"P„(y)) < 2-y"e. (4.50) 

Letting j/(e) be the largest value of y such that these inequalities simultaneously hold for infinitely many n, we 



conclude by Lemma 4.1C that 



y{e) < Y^^S{p\\a). (4.51) 
In particular, if y > S'(p||(t), then there exists e such that y > y(e), so 

limsupP„(y) < 1-e (4.52) 

n — >oo 

as required. □ 
Remark. Similarly, using the fact that Pn{y) is optimal among projections, we can conclude from the other half 



of Lemma 4.10 that lim„^oo Fn{y) = 1 when y < S{p\\a). We also have the natural conjecture that the lemma 
can be strengthened to say lim„^oo Fn{y) = when y > S{p\\a). 

This, of course, suggests that we should remove the requirement that a be p.p.t.; the same proof then gives: 
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Theorem 4.13. For any states p and a, 



Dr{p)<S{p\\a)+\og,Ti-\a^\. (4.53) 

When a is p.p.t., we recover the previous bound, while when a = p, we obtain the bound of |^ (see the 
remark following Corollary |4.6| above) . Note that we could also have obtained this result using Theorem 1 of 



pO| , based on the fidelity bound of Corollary 4^; this is essentially just the dual of the above proof. [| The proof 
given above was chosen to emphasize the fact that any bound on distillable entanglement can in principle be 
deduced from the dual SDP bound. 
If we define 

Bip,a) -.^ S{p\\a) + log^Trla^l, (4.54) 

then we have: 

Theorem 4.14. For any states p and a, and any trace-preserving superoperator ^ with both ^ and 4''" positive, 

B{^{p),^'{a))<B{p,cy). (4.55) 
For any other state p' and real number < p < 1, 

B{pp + (1 - p)p', a) < pBip, a) + (1 - p)B{p\ a). (4.56) 

Finally, we have in general 

B{p (dp', (J® a') = B{p, a) + B{p', a'). (4.57) 

Proof. Indeed, this is true for each of the functions S'(p||(t) and logjTrjcr'"! individually, so must be true for 
their sum. □ 

In general, B is not convex in a. In particular, we cannot assume that a local maximum of B is necessarily 
a global maximum. This is likely to make it very difficult to explicitly compute rmn^ {B{p, a)), although one 
can still, of course, obtain bounds from any given value of a. 



5 Exploiting symmetries 

If the state p has a large group of local symmetries, we can greatly simplify the primal and dual SDPs, in several 



cases to the point of being linear programs. The key observation is that, by the proof of Theorem 4.7, we have: 



Theorem 5.1. Let be a trace-preserving superoperator with both 4* and 4''" positive. Then for any state 
p — 4'(/o) and any K > 0, if F is primal optimal and D dual optimal for Fr{p;K), then so are 4'*(_F) and 
4'(Z?). In particular, if = 4*, we may assume that F is "^i* -invariant and D is ^ -invariant. 



^M., P., and R. Horodecki (personal communication) have pointed out a third proof via Theorem 2 of it is reasonably 
straightforward to show that the new bound satisfies their criteria for an upper bound to distillable entanglement. 
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Corollary 5.2. Let G be any closed subgroup ofU{k) <8) U{1), and let p be a G-invariant state; that is, a state 
such that for all U gG, 

UpU^ = p. (5.1) 
Then for any K > 0, there exists primal optimal F and dual optimal D invariant under G. If we further have 

Uop'Ul = p, (5.2) 



for some Uq e f/(fc) (8i U{1) with UqGUq = G, then we may further take 

UoF'ul = F, (5.3) 

UoD'U^o = D. (5.4) 

Proof. Let ^ be the superoperator 

^ :p^ I UpU\ (5.5) 

JU^G 

integrating with respect to the uniform probability measure on G. This is trace-preserving, e-local (thus p.p.t.), 
and satisfies vj/ = vf* = ^2 rpj^g ^^^^ claim thus follows from the theorem. 
Similarly, if is the superoperator 

^' ■.p^\{^{p) + Uonpful), (5.6) 

then the theorem applies to ^f'. □ 

Remark. In particular, if p is real, then we can take Uq = 1, allowing us to force F and D to be real as well. If 
p = Uop*^Uq for some Uo, we will say that p is pseudo-real. 

To apply this, it will be helpful to work in greater generality initially. Suppose simply that p is a Hermitian 
operator invariant under a subgroup G C U{k); we would like an efficient representation of p in which it is still 
straightforward to test positivity. 

Clearly, p is invariant under G if and only if p commutes with every element of G. But then p in fact 
commutes with the algebra C[G] of linear combinations of elements of G. In other words, p must be an element 
of the eentralizer algebra A of C[G']. From representation theory, we have: 

Lemma 5.3. There exists a unitary change of basis exhibiting an isomorphism 

C[G] ^ ®x (Mat(dA, C) ® 1„J , (5.7) 

for appropriate constants d\ and m\ such that 

^di = dim(C[G]), (5.8) 

A 

mxdx = k. (5.9) 

A 

In the same basis, the eentralizer algebra is given by 

eA(ld;, (8>Mat(mA,C)). (5.10) 
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In particular, the state p is determined by a set of Hermitian operators p\, with dimensions my, furthermore, 
p is positive if and only if p\ is positive for each A. Pseudo-reality conditions also carry over readily: in an 
appropriate basis, they produce conditions of the form (a) p\ real, (b) p\ = pv , or (c) p\ quaternionic. Finally, 
we have the trace identity 

TV(pa) = ^dATr(pA«TA). (5.11) 

A 

In particular, our simplification of Ft above can be viewed as a special case of this, based on the following 

two examples: 

Example. Let d be the subgroup of U{d^) consisting of operators U ®U. Any Gj-invariant operator can be 
written in the form 

p = x^{d)+y{l-^{d)), (5.12) 

with p > iff X, y > 0. 

Partial-transposing the above example, we get: 

Example. Let Gw be the subgroup of U{d'^) consisting of operators U ®U. Any G,(, -invariant operator can be 
written in the form 

p=|(l + rf$(d)r) + |(l-rf$(rf)r) (5.13) 

with p > iff a;, y > 0. 

Another important example is: 

Example. Let p be a state of dimension d. Then the state p*^" is invariant under the symmetric group Sn, acting 
by permuting the tensor factors. If p' is a generic 5'„-invariant operator, then the blocks p'^ are in one-to-one 
correspondence with the degree n representations of GLd(C), in such a way that p®" maps to the image of p in 
the corresponding representation. 

If p itself has symmetries, then we can simplify further. If A{x, y) is a homogeneous polynomial in two 
variables, then we write 

A{x,y))^Q (5.14) 
to denote the condition that A has nonnegative coefficients; similarly, 

A{x,y)<B{x,y) (5.15) 
means that B{x, y) — A{x, y) has nonnegative coefficients. 

Theorem 5.4. For any real numbers Q < f <1, K > Q and any integers d> 1, n> Q, 

Frihifr^; K) = maxi3(/, (1 - f)l{d^ - 1)), (5.16) 
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where B{x,y) and S{x,y) range over homogeneous polynomials of degree n such that 

O^B{x,y)^ix + {d^ ~l)yr (5.17) 

1 U^+d d^-dY ^, , I (d'^+d d^-dY /r.oN 

-x+^^y] ^S{x,y)^-{^—x + ^—y] (5.18) 



K \ 2 2 " - ^ - K 



5(.,,) = E(' (" + ^)-:(-^-^)^ ^L (5.19) 



Proof. Let F be primal optimal for Fr{Id{f)^^; K) such that F is invariant under Sn and G". The represen- 
tations of this group are in one-to-one correspondence with the integers < A < n, with d\ = (^) (d^ — 1)^ and 
m\ = 1. Writing 

Bx = dxFx, (5.20) 

B(a;,y)=^BAa;"-V, (5.21) 

we have 

Tr(F/rf(/)®") = d^P^^h (5.22) 

A 

= ^B,r-^((l-/)/(d2-l))^ (5.23) 

A 

= B(/,(l-/)/(rf^-l)). (5.24) 
We next observe that < iff ^ B{x, y) and F < 1 iff 

B{x, y)<Y. dxx'^S^ = {x + (rf2 - l)y)". (5.25) 

A 

Similarly, the partial transpose F^ is invariant under Sn and G^. Again the representations are indexed by 
< A < n, with 



Defining 



.Ay . 

= 1 (5.27) 



Sx = d'^iF^)x, (5.28) 
S{x,y) = Y,Sxx^-^y\ (5.29) 



we obtain the condition 



-x+^^y] ^S{x,y)^-{^—x+^—y] . (5.30) 



K \ 2 2 " - ^ - K 
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Finally, the relation between S{x,y) and B{x,y) obtains by noting that 



Six,y)^Tr[F^[^il + d^df) + la-dHdf)) ) (5.31) 
- Tr(F((l±l>i^i^a>(.) + ^(1 _ m)r) (5.32) 
^^.id+l)x-{d-l)y^x + yy (5.33) 

□ 

Remark. For d ~ 2, this linear program appeared in [Q, as an upper bound on the fidelity of separable 
distillation; the observation that it provides a lower bound on p.p.t. distillation is new. 

Similarly, 

Theorem 5.5. Fix a real number < p < I and an integer d > 2, and let Wd(j>) denote the Werner state 



Wdip) = ^^(1 + r(21)) + -^{1 - T(21)), (5.34) 



where T (21) = d^{d)^. Then 



whe 



F,(l^.(,)-;i.)=maxi.(f^,-^), (5.35) 



, , fd^ + d d^-d Y , , 

0:<B{x,y):< i-^x+^—y\ (5.36) 

^^{x + id' - l)yr ^ Six, y) ^ -i(x + id' - I)?/)" (5.37) 

Bix,y)^s( ^'^'^^-^'-'^\ ^-±l]. (5.38) 



Corollary 5.6. For the antisymmetric Werner state IFd(l), we have 



Fr{Wdil);K) = min(l, ^). (5.39) 
Dr{Wa{l)) = log2 . (5.40) 



For any state p, 



dK 

Frip®Wdiiy,K)^Frip;j-^). (5.41) 
Drip ^ Wdil)) = Drip) + DriWdil)). (5.42) 

Proof. We observe that Tr |W<j(l)^| = Thus if we show that Fr(W^d(l); ^) = 1, the proof of Corollary 

4.3 will apply to give equation ( |5.4l| ); taking p — <f>(l) gives equation ( 5.39 ), and the equations for Dr follow 



immediately. It thus remains to show FriWdil); ^^) ^ 1 (since the other inequality is immediate). 
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Taking 



B{x,y) 



S{x,y) 



d-2 
d + 2 



(P -d 



-y, 



d + 2 



{~x+id^ -l)y), 



we find 



Fr{Wd{l);^)>B{0, ^ 



d^-d 



= 1. 



(5.43) 
(5.44) 

(5.45) 
□ 



Similar results apply to "iso- Werner" states — states which are linear combinations of 1, T(21), and $(c?) 
(invariant under O^O with O G 0(d)) — and Bell-diagonal states — states on C^^^ which are linear combinations 
of $(2) and CTu,<i>(2)cr^^ for w G {x,y,z}. 

Using Theorem 4.14 , we can apply the argument of Corollary 5.2 to conclude that when minimizing B{p, a), 
we may insist that a possess the symmetries of p. When p is isotropic, we learn nothing new (the earlier bound 
(01 fl^li 0) unchanged), but when p is Werner, we obtain: 



Corollary 5.7. Fix a real number < p < I and an integer d > 2. Then 



DriWdip)) < mmB{Wd{p),c7) = I I + p\og^{p) + (1 -p)log^{l-p) i <p < i + i 



(5.46) 



\og,{^)+p\og, (^) 
Proof. By the above argument, we may assume a — Wd{p')- Now, 

B{Wdip), Wdip')) = plog, +{1-P) log, (^|-^ 

We find that the optimal p' satisfies 




(5.47) 



P 



P 

1 

2 



<p < 5 
h<P<h 



(5.48) 



^. d+2-4p 2 ^ d — — 

Plugging in, we obtain the stated bound. □ 

Remark 1. We observe that this bound is differentiable and convex for < p < 1, and tight for p = 1. 

Remark 2. The above bound has recently been independently derived in , as the regularized relative entropy 
of entanglement; that is, 



lim - min 5(p®"||(7). 

n^cso Jl (7 p.p.t 



(5.49) 



This suggests that the bounds of Theorems 4.11 and 4.13 may regularize to the same bound. 
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6 Hashing analogues 

One of the few known lower bounds on distillable entanglement is based on the "hashing" protocol |^ ; it will 
be instructive to consider this bound (for p.p.t. distillation) via the present techniques. The key point of the 
hashing bound is that on "low weight" states, it gives fidelity close to 1, while on "high weight" states, it gives 
fidelity close to 0. This suggests the reasoning behind the following proof: 

Theorem 6.1. Fix a fidelity | < / < 1 and an integer d > 1. Then 

i?r(/d(/)) >max(log2d + /log2/ + (l-/)log2i-4,0). (6.1) 

a+1 

Proof. Fix an integer n > 0, and consider the set Vn consisting of tensor products 

P = ®l<^<nP^ (6.2) 

with each Pi G {$((i), 1 — $((i)}; note, in particular, that Vn is a set of mutually orthogonal projections. Since 

Mdf\-\ and l^<|l-ci,(d)r|<l±l, (6.3) 

we have 

|P^| < 1)^*(^\ (6.4) 

where we define wt(P) to be the number of factors equal to 1 — <i>(rf). 
Let us then define an operator 



Wt(P)<lO 



We observe that Fn{w) is a projection, so < Fn{w) < 1, and that 

Tr(F„(z.)(/,(/)®") = ^ f")r-Xl-/r, (6.6) 



0<i<w ^ ' 

which tends to 1 as n ^ cx) as long as 

uj := lim wjn > 1 — /. (6.7) 

n — ^00 

We also compute 



PePn 0<i<W 
Wt(P)<U) 



Fn{wf\< J2 \P^\<d-'' E [Md+iy- (6.8) 



If we take uj < ^q^, then we obtain the limit 

lim ilog2(d-" V (''){d+iy)^-ujlog2Uj-{l-Lo)\og2{l-uj)+ujlog^{d+l)-\og^{d). (6.9) 



0<j<iu 
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But then by Theorem 3.1, we conclude that 

^r(/d(/)) >log2(d)+c^log2 



+ (1 - w)log2(l -W) 



(6.10) 



whenever 1 — f < uj < Since this is decreasing over the range, we obtain the strongest bound by taking 



the hmit as lj ^ 1 — /, proving the theorem. 



□ 



Remark. When d — 2, this is precisely the hashing lower bound (albeit weaker, in that it applies only to p.p.t. 
distillation). However, for d > 2, the new bound is strictly stronger. 

This gives us a general technique for proving lower bounds on p.p.t.-distillable entanglement: approximate 
the given state as a linear combination of projections with well-controlled partial transposes. Our primary 
application of this will be to "maximally correlated states" [|lOj . We recall that a maximally correlated operator 
is one of the form 



l<i,j<*: 



(6.11) 



for some positive Hermitian operator a, and similarly for a maximally correlated state. In [|T^, an upper bound 
was given on the p.p.t. distillable entanglement, and the conjecture was made that this bound was tight (even 
for the 1-locally distillable entanglement). We give a partial resolution of this conjecture: 



Theorem 6.2. For any maximally correlated state 

Pa = ^a.y|M)(jj|, 



the p.p.t. distillable entanglement is given by the formula 

DriPa) = B{a) := H{aii,a22, ■ ■ ■) - S{a) 



(6.12) 



(6.13) 



Proof. That this is an upper bound was shown in , so it suffices to prove the lower bound. We construct a 
protocol in two steps. 

First, suppose a possesses a transitive group of symmetries; that is, a transitive group G of permutations 
such that 



(For instance, the operator 



a7r(i)7rO) = a^j, Vtt G G, 1 <i,j < k. 



6 1 



(6.14) 



/l 

3 


i 
6 


— i 


1 


6 


3 


u 


— i 
6 



(6.15) 



is symmetric under the transitive group Z3 of cyclic shifts.) We decompose 

^ ApA(n), 

AG(0,1] 



(6.16) 
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where px{n) is the orthogonal projection onto the A-eigenspace of a®". Then p\{n) is symmetric under the 
transitive group G", and thus has constant diagonal. If we similarly decompose 

pr= E ^^^N' (6.17) 

Ae(o,i] 

we find 

Px{n) = Pp^(n)- (6.18) 
We can thus apply the following lemma to Px{n). 

Lemma 6.3. Let p = be a maximally correlated operator of dimension d x d such that (3 has constant 
diagonal. Then 

l/l < (6.19) 

Proof. We compute 

p'^ = EAi|ii>(ii|- (6.20) 

This is a block matrix with 1- and 2-dimensional blocks: wc thus immediately compute that its eigenvalues are 
Pa ioT 1 < i < d, and ±|/3y| for 1 < i < j < d. Since /3 is positive, we have: 

<PuPn, (6.21) 



and thus the largest eigenvalue of p^ in absolute value is maxi<i<ci Pa. We thus have 

_ Tr(p) 

l<i<d 



max I3u = l V 0^^ = (6.22) 
l<i<d d ^ d 



□ 



Now, write 

P{n,x)= E ^A(n), (6.23) 

x<A<l 

Xe = mi{{x G (0, 1] : Tr{pP{n, a;)) > 1 - e}). (6.24) 



Then for e > 0, we find that since 



we have 



Drip) > Hm hm -- log^ ^(Pi^^^^)) (g 26) 

£^0n^oo n a" 

= \og^d-S{a). (6.27) 



22 



Since H(aii, a^i . . . ) = H{\/ l/d, . . . ) = logj d, we have proved the theorem in the symmetric case. 

To reduce the general case to the symmetric case, we adapt the distiUation protocol for pure states given in 
[Q. Given a word w in the numbers 1 ... fc, we write wti(u') for the number of times i appears in w. Then our 
first step is, given 

pT=Y.\ n I (6-28) 



W.W 



, l<m<r; 



to measure wt? for 1 < i < k. Then the resulting (random) state p^' is maximally correlated, and a' admits a 
transitive action of 5'^. Now 

Dripa) = -DripT) > - E(i?(p„0) - - E(S(a')), (6.29) 
n n n 

where E(-) is the expected value, and the inequality follows from the fact that the measurement is local, so 
cannot increase the expected distillable entanglement. It thus suffices to show that 

E(S(a')) = nB{a) + o{n). (6.30) 

Now, the measurement has at most n'^ different outcomes, so gives us at most k log2 n bits of information. But 
then 

E{H{a')) > nH{a) - k log^ n, (6.31) 

E(S'(a')) < nS{a), (6.32) 

so we find 

E{B{a')) > nB{a) - klog^n = nB{a) + o(n), (6.33) 

as required. □ 

We also have the following general result. 
Theorem 6.4. Fix a finite- dimensional Hilbert space V , and let 

l<i<m 

he a partition of the identity with the Pi orthogonal projections. For each 1 < i < m, let Si be the largest 
eigenvalue of \P[\. Then for any state p £ V{V ® V), we have 

Dr{p)> ^ r,;(log2r, -log2Sj), (6.35) 

l<i<m 

where 

r, := Tr(P,;p). (6.36) 
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Proof. To any subset 5 C {1, 2, . . . m}", we associate a projection 

Ps^J2 (6.37) 

toes l<'i<n 

which satisfies 

\Ps\<T. n (6.38) 

w£S l<i<n 

For each < e < 1 and each integer n > 1, let /3„(e) be the minimum over S of the largest eigenvalue of |Pg | 
subject to the constraint Tr(Psp'^") > 1 - e. Then 

Drip) > liminfliminf ^^^^^^^ (6.39) 

(take F ^ Ps). Since 

Tr(P5P^") = E n (6-40) 

lueS l<i<n 



the theorem follows by the classical analogue of Lemma |4T0|. □ 



7 Clones 

In this section, we sketch a possible direction to take in applying the above techniques to 1-local questions 
(quantum codes and distillation protocols). 

Definition 5. An operator A on (C'')®" is an "n-clone" if it can be written in the form 

A = Y.Af- (7.1) 

i 

where each Ai is a positive operator, or can be written as a limit of such operators. 

For a permutation tt e S'„, T{n) is the operator on (C*^)®" that permutes the tensor factors by tt; when 
TT = 21 G 52, this agrees with our earlier notation. 

Theorem 7.1. Let A be an n-clone. Then for all involutions n g Sn, and all sets S C {1, 2, . . . n} that intersect 
each 2-cycle of t: exactly once, the following operator is positive: 

{ATin)f^. (7.2) 

Proof. Since nonnegative linear combinations and limits of positive operators are positive, it suffices to prove 
the result for A — A®". In that case, {AT{Tr))^^ factors as a tensor product of the following operators: 

A, A\ and {{A (g)A)T (21) f\ (7.3) 

The first two are clearly positive; that the third is positive is a special case of the following lemma. □ 
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Lemma 7.2. For any operator A (not necessarily Hermitian) , the operator 

is positive. 
Proof. We have 

{{A ® A^)T{2l)f- - {{A ® l)r(21)(At ® 1))^^ ^ (A ® l){T{2l)f\A^ ® 1). 
Since T{2lY^ = dim(y4)$(dim(A)) > 0, the resuh follows. 



(7.4) 



(7.5) 
□ 



For instance, let C be a quantum code of length n over an alphabet of size k, and consider the following 
average over codes equivalent to C: 



W{C) = ^c'^c{Pc' ®Pc'). 
This is clearly a 2-clone, so we conclude that the following operators are positive: 

W{C), W{Cf, {W{C)T{2\)f. 



(7.6) 



(7.7) 



We also find that W{C) is invariant under operators of the form U ®U , with U in the semidirect product of 
Sn acting on U{k)^'^. Thus using the techniques of Section ||, we conclude that the three given operators are 
positive if and only if the following three polynomials have nonnegative coefficients: 



Sc{x,y) ■.^TriW{C)ix 



1-T(21) l + r(21). 



Bc{x, y) :- Tv{W{CY {{x-T{2l) + y(l - -r(21)))«")^ ) 



Ac{x, y) := Tv{{W{C)T {2l)Y {{x-T{2l) + y{l - -r(21)))«")^ ). 



Using the fact that Ty{M^N^) = Tr(MiV), we find 



where 



5c(x,y)=A^(^,^), 

Bc{x,y) = ^ciy^^). 
Acix,y)^A'c{^^,y), 



A'cix.y) :=Tr(Ty(C)(.T + yT(21))«"). 



(7.8) 
(7.9) 
(7.10) 

(7.11) 
(7.12) 
(7.13) 

(7.14) 



In other words, these are precisely the weight enumerators of C ([p^, pT|). In the full linear programming 
bound for quantum codes, there is an additional inequality: 

1 



Bc{x,y) 



dim(C) 



Ac{x,y) > 0. 



(7.15) 
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To prove this, we simply extend C to a self-dual code C+ by encoding half of <i>(dim(C)) into C . We then have 

A'c+{x,y,u,v) = A'(.{x,y)u + A'(j{y,x)v, (7.16) 



c / N S{x,y) ~ S{~x,y) S{x,y) + S{-x,y) 

Sc+ [x, y, u, v) = u H V (7.17) 

Bc+ {x, y, u, v) = ^.^^^^^ ^c(a^, y)u + {Bc{x, y) - ^.^^^^^ Ac{x, y))v (7.18) 

Ac+ [x, y, u, v) = '^^~((j^^c{x, y)u + {Bc{x, y) - ^.^^^^^ Ac{x, y))v. (7.19) 

In particular, the polynomial Bc{x,y) — -^^r^^j^ Ac {x , y) must have nonnegative coefficients. 

We can thus extend the linear programming bound to higher-order invariants (p^) by using the relevant 
symmetry group to decompose the operators attached to 

by Theorem Note that since Wi (C+)r(7r) — Wi (C+) for tt e S'„, we have only [|J + 1 operators to consider. 

Another application of the clone concept is to 1-local operations. Fix a Hilbert space Va<8Vb, and consider 
the 1-local operation 

^ = ^A®B,, (7.21) 

i 

where Bi are operations, and Ai are completely positive superoperators such that ^ - Ai is an operation. Then 
as remarked in H] , we can extend to an operation on the larger Hilbert space Va ® Vg"' by simply taking 

^(") =^A,«)Sf". (7.22) 



Note that this depends not just on ^E* but also on the specific decomposition (7.21). The following is straight- 
forward: 

Lemma 7.3. For any 1-local operation ^! , any integer n > 1, and any vector u € Vyi Va, the operator 

TvAmivl ® l)(*("H<i>(l^A ® Vi"m (7.23) 

is an n-clone. 

Using Theorem 7.1, we obtain a number of semidefiniteness constraints that must satisfy; these con- 
straints can in principle be used to obtain bounds on 1-local distillation. (For instance, the argument of ^ can 
be restated in these terms, although we have not done so.) Unfortunately, the resulting semidefinite programs 
tend to be fairly complicated, and thus further ideas would seem to be needed. We also note that the cloning 
argument is quite fragile; if we define a "catalyzed" fidelity 

Fi{p-K) =limsupFi(p0$(d);dii'), (7.24) 

d—^OG 
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after Corollary |4.3| , then we can no longer directly use cloning to bound the corresponding distillable entangle- 
ment. 

We close with the following new application of the cloning argument: 

Theorem 7.4. Fix a pair of integers I < K < d. Then for all fidelities 1/d < f < 1, we have the strict 
inequality 

FMf);K) < 1/K + l/K). (7.25) 

Proof. Suppose we had equality. A protocol ^' attaining this bound would certainly have to be p.p.t.; thus, if 
we apply this protocol to Id{f'), the output fidelity will take the form F{f') ~ af + b for some constants a and 
6, or equivalently 

for constants a', b' . Evaluating this at /' = l/d, /' = 1, we find: 

a' <l/K, b' <l. (7.27) 

On the other hand, at /' = /, we have 

(7.8, 

Since the coefficients are both positive, we conclude that a' — l/if, b' = 1. In particular, 5* must take 
Id{l) - $(rf) to Ik{1) - ^{K). 

Now, consider the action of ^P^^) on the state $(d)® ^1^. Since ^ takes the pure state <l>(d) to the pure state 
(f>(if ), we conclude that must take ^(c?) ® ^1^ to a state of the form ^{K) ®X] by symmetry, we conclude 
that X — -^Ik- But then tracing away the other copy of Vb, we find that ^E" takes Id{l/d^) to Ik{1/K^)- On 
the other hand, we have 

We thus obtain a contradiction, and the theorem follows. □ 
Remark. From ||^, it follows that 

FMf);K) > l/K' + ^^-^{l - l/K') (7.30) 
whenever < K < d. Is this lower bound tight? 
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